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Abstract 

Bats may be natural reservoirs for a large variety of emerging viruses, including mammalian coronaviruses (CoV). The 
recent emergence of severe acute respiratory syndrome-associated coronavirus (SARS-CoV) and Middle East respiratory 
syndrome coronavirus (MERS-CoV) in humans, with evidence that these viruses may have their ancestry in bats, highlights 
the importance of virus surveillance in bat populations. Here, we report the identification and molecular characterization 
of a bat B-Coronavirus, detected during a viral survey carried out on different bat species in the island of Sardinia (Italy). 
Cutaneous, oral swabs, and faecal samples were collected from 46 bats, belonging to 15 different species, and tested for 
viral presence. Coronavirus RNA was detected in faecal samples from three different species: the greater horseshoe bat 
(Rhinolophus ferrumequinum), the brown long-eared bat (Plecotus auritus), and the European free-tailed bat (Tadarida 
teniotis). Phylogenetic analyses based on RNA-dependent RNA polymerase (RdRp) sequences assigned the detected CoV 
to clade 2b within betacoronaviruses, clustering with SARS-like bat CoVs previously reported. These findings point to the 
need for continued surveillance of bat CoV circulating in Sardinian bats, and extend the current knowledge on CoV ecology 
with novel sequences detected in bat species not previously described as B-Coronavirus hosts. 


Keywords Bats - Coronavirus - RNA-dependent RNA polymerase - Rhinolophus ferrumequinum - Sardinia 


Introduction high diversity and ecological characteristics that potentially 
facilitate virus transmission, such as the roosting together 
behavior and their ability to fly over distance [4, 5]. There- 


fore, it is essential to identify and characterize the currently 


Bats (order Chiroptera) are one of the most diverse and 
widely distributed group of mammals, comprising more 


than 1300 species across 21 families [1], and play a key 
role in terrestrial ecosystems [2]. In the last few decades, 
they have been associated with several emergent zoonotic 
agents and recognized as potential reservoir of zoonotic 
viruses belonging to different families, such as Hendra, 
Nipah, Ebola, Lyssa, and SARS Coronavirus [3]. Bats rep- 
resent ideal hosts and reservoirs of viruses, due to their 
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unrecognized viruses circulating in bat populations. 

Coronaviruses (CoV, family Coronaviridae) are envel- 
oped viruses with large single-stranded RNA genomes rang- 
ing between 27 and 32 kb [6, 7]. Based on genome-scale 
phylogenies, they are classified into four genera: Alphacoro- 
navirus (a-CoV) and Betacoronavirus (B-CoV), pathogenic 
for mammals, Gammacoronavirus (y-CoV) and Deltacoro- 
navirus (6-CoV), primarily (but not exclusively) circulating 
in birds [8, 9]. Most CoVs are disease-causing agents, asso- 
ciated with a range of respiratory, gastrointestinal, hepatic 
and neurological diseases both in humans and animals. They 
are also potentially transmitted among species, causing epi- 
demics after host transition [10—12]. Due to mutation and 
recombination events, new strains constantly originate and 
spread [4, 13]. 

Different bat species have been identified as natural reser- 
voirs of most mammalian a-CoV and B-CoV, including the 
Severe Acute Respiratory Syndrome SARS-CoV, and the 
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Middle East Respiratory Syndrome MERS-CoV [14, 15], 
which recently caused outbreaks in the human population 
[4]. SARS-related CoV was first identified in Rhinolophus 
bats from China [16], and subsequently many mammalian 
CoVs were detected in Chinese bats [17, 18] and globally 
[19]. Further investigations identified a-CoV and B-CoV 
in several European countries from bats of various species 
[20-22]. 

In Italy, few studies have been conducted on bat CoV 
ecology: viral detection and prevalence was investigated in 
greater horseshoe bats [23], alpha- and betacoronaviruses 
were identified in different species, such as Kuhl ’s pipist- 
relle, lesser horseshoe bats and serotine bats [10, 24], and 
recently two MERS related beta CoV detected in Italian 
Kuhl’s pipistrelle and Savi’s bat were fully sequenced [25]. 
A total of 21 bat species are currently recognized as liv- 
ing on the island of Sardinia (Italy), lying geographically 
at the interface between Europe and Africa, and recently 
indicated as a major Mediterranean hotspot for bat biodi- 
versity [26]. Nothing is known about CoV presence and 
distribution among bat populations on the island, and this 
is the first study aiming to fill this gap. Herein we report the 
molecular identification and typing of a SARS-like Corona- 
virus detected in three species of Sardinian bats belonging to 
Rhinolophus, Plecotus and Tadarida genera, during a viral 
survey conducted on non-invasive bat samples collected 
from different localities on the island. We also inferred the 
phylogenetic characteristics of the identified B-CoV, and its 
evolutionary relationships with other bat CoVs identified in 
Europe and globally. 


Table 1 List of bat samples 
analyzed: host species, sampling 
locations, and type of samples 
collected 


Host species 


Rhinolophus mehelyi 
Rhinolophus ferrumequinum 
Miniopterus schreibersii 
Myotis punicus 

Myotis capaccinit 
Myotis emarginatus 
Plecotus sardus 
Hypsugo savii 
Pipistrellus pipistrellus 
Pipistrellus kuhlii 
Plecotus auritus 
Plecotus austriacus 
Nyctalus leisleri 
Eptesicus serotinus 
Tadarida teniotis 

Tot 
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Materials and methods 
Bat sample collection 


A total of 46 bats, belonging to 15 different species, were 
captured using mist nets or hand nets in 8 different locali- 
ties distributed in northern and central Sardinia: Grotta sa 
Rocca Ulari (Borutta), Galleria Casteldoria (Santa Maria 
Coghinas), Diga Santa Chiara (Ula-Tirso), Monte Coral- 
linu (Dorgali), Tilipera (Bonorva), Mularza Noa (Bolo- 
tana), Belvi, and Grotta Verde (Alghero), listed in Table 1. 
Sampling locations included caves, buildings, tunnels, 
woods, abandoned bridges, and other type of natural roosts. 
Captures were performed by expert and trained biologists 
within a regular census program of Sardinian bat popula- 
tions, between November 2015 and November 2016. Bats 
were placed individually into cotton bags before investiga- 
tion. Species, sex, age category, and forearm length were 
determined; oral and cutaneous swabs were collected using 
sterile micro-swabs (Aptaca), and then bats were released 
at their capture site. Faecal samples were readily collected 
from cotton bags or during bat handling and suspended in 
1 ml RNA later RNA stabilization solution (Qiagen) and 
then stored at — 20 °C before processing. 


DNA and RNA extraction, reverse transcription 


DNA was extracted from oral and skin swabs tips using 
the QiAmp DNA Mini Kit (Qiagen), following the 


Sampling location Skin swab Oral swab Faeces 
1 2 — 2 
1,2,8 3 3 3 
1 9 — T 
1,2,5 8 8 6 
2,8 3 4 3 
2,6 3 2 2 
3 — 1 — 
4,6 1 1 — 
5 ps — — 
6 — — 1 
7 1 1 1 
8 — 1 — 
6 1 1 — 
8 — — 1 
8 1 1 1 
34 23 27 


Locations: (1) Grotta sa Rocca Ulari (Borutta), (2) Galleria Casteldoria (S.M. Coghinas), (3) Diga S.Chiara 
(Ula-Tirso), (4) Monte Corallinu (Dorgali), (5) Tilipera (Bonorva), (6) Mularza Noa (Bolotana), (7) Belvì, 


(8) Grotta Verde (Alghero) 
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manufacturer’s protocol specific for swabs. DNA was stored 
at — 20 °C for subsequent amplification analysis. After sus- 
pending faecal samples in phosphate buffered saline (PBS), 
RNA was extracted from faeces by using the QiAmp Viral 
RNA Mini Kit (Qiagen), according to the instructions sug- 
gested by the manufacturer. Extracted RNA was stored at 
— 80 °C. First-strand cDNA was synthesized utilizing the 
Affinity Script Multi Temperature cDNA Synthesis Kit (Agi- 
lent Technologies), following the random primers strategy 
protocol, and then stored at — 20 °C. 


PCR and RT-PCR, screening and sequencing 


DNA extracted from bat oral swabs was tested for the pres- 
ence of Herpesvirus using a universal nested PCR [27]. 
DNA extracted from skin swabs was screened for Poxvirus 
and Papillomavirus presence, by PCR based on pan-pox 
primers [28] and on primers targeting a conserved region 
of PPV LI gene [29]. 

RT-PCR was performed to detect Coronavirus in bat 
faecal samples, by using two pairs of primers targeting a 
conserved fragment of RNA-dependent RNA polymerase 
gene (RdRp): 11FW and 13RV [23], and Corona | forward 
and Corona 2 reverse [10, 17, 30], targeting 170 bp and 
440 bp respectively. In a typical PCR, 2 ul of cDNA was 
amplified in a 50 ul reaction mixture containing 0.2 mM 
of each dNTPs, 0.2 uM each of Forward and Reverse prim- 
ers, 1X TaqBuffer and GoTaq (Promega, Italy) DNA poly- 
merase. The cycling conditions were as follows: 5 min at 
95 °C followed by 40 cycles of 1 min at 95 °C, 1 min at 
50 °C, and 1 min at 72 °C, followed by a final extension 
time of 5 min at 72 °C. From the first amplified sequence 
of approximately 170 bp, in order to further amplify the 
RdRp gene, 4 sets of primers were newly designed, after 
aligning the sequence detected in this study with the most 
similar sequences obtained by BLAST analysis, and choos- 
ing conserved regions to design consensus walking prim- 
ers. The first 170 bp amplicon was extended by design- 
ing and using downstream primers within the fragment 


Table 2 Primer sequences used 


o Name Target sequence 

and designed in this study for 

bat Coronavirus detection 11FW RdRp 
13RV 
BatCV1F RdRp 
BatCVIR 
BatCv2F RdRp 
BatCV2R 
BatCV3F RdRp 
BatCV3R 
BatCV4F RdRp 
BatCV4R 


and upstream primers based on alignments of phyloge- 
netically related CoV from GenBank. Table 2 lists the 
primers designed in this study (BatCV1F/R, BatCV2F/R, 
BatCV3F/R, and BatCV4F/R). PCR products were purified 
using DNA Clean & Concentrator purification kit (Zymo 
Research, CA, USA), and directly sequenced in both for- 
ward and reverse directions by BMR Genomics Sequenc- 
ing Service (Padova, Italy). Forward and reverse sequences 
were edited and aligned with MUSCLE Multiple Sequence 
Alignment tool (http://www.ebi.ac.uk/Tools/msa/), align- 
ments were edited with BioEdit software [31]. Sequences 
were checked against the GenBank database using nucleo- 
tide BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cg1). 


Phylogenetic analysis 


Phylogenetic trees were estimated using the Maximum 
Likelihood (ML) and Maximum Parsimony (MP) meth- 
ods implemented in MEGA 6 [32], after testing the best 
fitting evolutionary model to the data. Bootstrap support 
values and reliability of internal branches were calculated 
from 1000 replicate trees. Models with the lowest BIC 
scores (Bayesian Information Criterion) are considered to 
describe the substitution pattern the best. Non-uniformity 
of evolutionary rates among sites was modeled by using 
a discrete Gamma distribution (+G) with 5 rate catego- 
ries and by assuming that a certain fraction of sites are 
evolutionary invariable (+I). After generating a gap-free 
nucleotide alignment of 1090 bp containing the novel virus 
as well as of Coronavirus RdRp reference strains, the most 
likely nucleotide substitution model selected to analyze 
these data was TN93+G (Tamura-Nei, Gamma distrib- 
uted), while analyzing an alignment of a higher number 
of CoV RdRp sequences of 396 bp, the model T92+G was 
selected (Tamura 3-parameter, Gamma distributed). The 
CoV sequence obtained in this study was submitted to the 
GenBank database under accession number MG978754. 


Primer sequence Source 


5'-TGATGATGCCGTCGTGTGCTACAA-3' 
5'-TGTGAGCAAAATTCGTGAGGTCC-3’ 


Balboni et al. [23] 


5'-TACAGTGATGTAGAAACTCC-3' This study 
5’-AAGGTCTGTCTCAGTCCAGCA-3' 
5'-TGCTGGACTGAGACAGACCTT-3’ This study 
5'-GTGTTATCATTAGTTAGCATT-3’ 
5'-CTAAGCGCAATGTCATCCCT-3' This study 
5’-GCATGTTAGGCATGGCTCT-3’ 
5'-GCATCAAGAACTTTAAAGCAG-3’ This study 


5’-GAATACATGTCTAACATGTG-3' 
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Results 


Coronavirus RNA was detected in three out of 27 faecal 
samples: 11%, close to previously reported frequency [10], 
less than that observed in other studies [17, 30], unlike 
results obtained in central-southern Italy [33]. Positive 
samples were collected in different locations of northern 
and central Sardinia from three different bat species: the 
greater horseshoe bat (Rhinolophus ferrumequinum), the 
brown long-eared bat (Plecotus auritus), and the European 
free-tailed bat (Tadarida teniotis). Table 3 lists samples 
positive for CoV detected in this study, including the Gen- 
Bank accession number assigned to the novel sequence. No 
Herpesvirus, Poxvirus and Papillomavirus were detected 
in the oral and/or cutaneous samples collected within this 
Sardinian bats viral survey. 

After detecting CoV in faecal samples, a total of 
1110 bp of RdRp gene were successfully amplified and 
directly sequenced. The same betacoronavirus sequence 
was detected in all three positive bat species. Upon 
BLASTN analysis, the viral sequence identified was found 
highly similar to Bat coronavirus BM48-31/BGR/2008 
(GU190215), detected in a Blasius’s horseshoe bat (Rhi- 
nolophus blasii) in Bulgaria [14], with 100% query cover- 
age and 96% identity. The obtained sequence also showed 
high nucleotide similarity (89%) with Bat SARS CoV 
Rs672/2006 (FJ588686) identified in a Chinese rufous 
horseshoe bat (Rhinolophus sinicus) [34], and with Bat 
SARS-like coronavirus RSSHCO14 (KC881005) from the 
same bat species [35]. A first phylogenetic analysis was 
performed based on 1090 bp amplified fragment of the 
RdRp gene including the CoV identified in Sardinian bats 
and other 16 sequences representatives of CoV genetic var- 
iability. The phylogram (TN93+G, Fig. 1) clearly showed 
that Sardinian bat CoV belong to the B-CoV genus, and 
clustered together with the Bulgarian bat CoV identified in 
a Rhinolophus blasii (GU190215). To have a deeper look 
into clade structuring within genera and further character- 
ize the phylogenetic relationships of the detected CoV, 
especially with other European bat coronaviruses, another 
phylogenetic analysis was performed based on 396 bp of 
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the RdRp gene, and including a total of 30 sequences. Fig- 
ure 2 shows the (T92+G) tree obtained, clearly assigning 
the Sardinian bat CoV to clade 2b within the B-CoV genus. 
The branching also showed that the detected sequence was 
closely related to a betacoronavirus previously identified 
in a greater horseshoe bat (Rhinolophus ferrumequinum) 
in Italy (BtCoV/Rhi_fer/It17/ITA/2009, GenBank number 
KC633199), as well as other bat CoVs detected in Rhi- 
nolophus spp. from Slovenia (GQ409747) [36], and again 
from Italy G.e., KF500954) [10]. All other clade groupings 
were consistent with previous studies [14]. 


Discussion 


At least 21 bat species are currently recognized as living 
on the island of Sardinia, Italy [37]. Some species are com- 
mon and adaptable, inhabiting both wild and urban habitats, 
while others are rare and restricted to limited areas of the 
island. Recently an island endemic species was character- 
ized (Plecotus sardus) and an endemic lineage of Plecotus 
auritus was also identified [38]. The existence of endemic 
bat taxa on Sardinia highlights the island importance in the 
conservation of the European bat community. Some Sardin- 
ian bat populations appear genetically differentiated and iso- 
lated from any continental gene pool and should therefore be 
considered as an evolutionarily significant unit (ESU) [39]. 

This is the first report of Coronavirus detection in faeces 
from bat species in Sardinia. A SARS-like betacoronavirus 
was detected in three species (Table 3): the greater horse- 
shoe bat (Rhinolophus ferrumequinum), common cave spe- 
cies distributed all over the island; the brown long-eared 
bat (Plecotus auritus), rare species inhabiting few locali- 
ties (especially mountain woods) in central Sardinia; and 
the European free-tailed bat (Tadarida teniotis), common 
species usually colonizing rocky slopes and Mediterranean 
bushes as well as buildings and bridges in urban areas. 

The detected virus (SarBatCoV1) belonged to clade 2b, 
including most of the SARS-like CoVs found in bats, as 
well as human SARS Coronavirus (Figs. 1, 2). SarBatCoV1 
appeared related to another betacoronavirus previously iden- 
tified in the Italian greater horseshoe bat (BtCoV/Rhi_fer/ 


Table 3 List of bat positive samples for Coronavirus, and GenBank accession number 


ID Host species Sample type 
15F Greater Horseshoe bat R. ferrumequinum Faeces 
30F Brown Long-eared bat P. auritus Faeces 
31F European Free-tailed bat T. teniotis Faeces 


Sex Location CoV strain GenBank 
accession 
number 

F 2 SarBatCoV 1 MG978754 

nd 7 SarBatCoV1 — 

F 8 SarBatCoV 1 — 


Locations (ref. Table 1): 2. Galleria Casteldoria (S.M. Coghinas), 7. Belvì, 8. Grotta Verde (Alghero) 


A Springer 


Virus Genes 


Bat SARS CoV Rs672/2006/R.sinicus/F J588686 
Bat SARS coronavirus Rp3/2004/DQ071615 
SARS-related bat CoV/R.sinicus/KF294440 
Coronavirus Lyra1 1/R.affinis/KF569996 
100 | Bat SARS-like CoV RSSHC014/R.sinicus/KC881005 
BtRs-BetaCoV/YN2013/R.sinicus/KJ473816 
29 | | SARS-like CoV WIV16/R.sinicus/KT444582 
Bat SARS coronavirus Rm1/R.macrotis/DQ412043 Beta 

100 Bat CoV isolate Jiyuan-84/R.ferrumequinum/KY 770860 


69 
100 “— Bat SARS coronavirus Rf1/R.ferrumequinum/DQ412042 


60 


92 


51 


50 SarBatCoV1 
87 — Bat CoV BM48-31/BGR/2008/R. blasii/GU190215 
Bat coronavirus HKU9-3/EF065515 
100 Bat coronavirus HKU9-4/EF065516 
Bat coronavirus HKU4-1/EF065505 


Bat coronavirus HKU6/DQ249224 | 
Alpha 


100 Bat coronavirus BtCoV/512/2005/DQ648858 


0.1 


Fig. 1 Phylogenetic analysis based on 1090 bp of the RdRp gene. 
RNA-dependent RNA polymerase-based phylogeny including the 
novel bat coronavirus detected in this study (SarBatCoV1) and other 
16 RdRp sequences from bat coronaviruses. Coronavirus sequence 


It17/ITA/2009, GenBank number KC633199), and also 
clustered with a B-Coronavirus isolated from a Blasius’s 
horseshoe bat in Bulgaria (BM48-3 1/BGR/2008, GenBank 
number GU190215). As shown in Fig. 2, the two clades 
within B-Coronaviruses are strongly supported in the phylo- 
gram, and clade 2c clustered many other bat CoVs identified 
in Europe (Italy and Spain) from different species. These 
findings confirm a correlation between SARS-like CoVs and 
bats belonging to Rhinolophus spp., consistently with previ- 
ous studies [14, 40], although the same RdRp sequence, in 
this study, was also detected in two other species belong- 
ing to the genera Plecotus and Tadarida. As previously dis- 
cussed, phylogenetic analyses providing insights into viral 
transmission among bat species have revealed a general clus- 
tering by geographic location rather than by bat species [4]. 
However, our results also show that similar coronaviruses 
(with nucleotide identity up to 96%) are found in bats dis- 
tributed in geographically distant regions, such as the island 
of Sardinia, mainland Italy and Bulgaria. 

Since the identification of SARSr-Rh-BatCoV in bats 
of the Rhinolophus genus [16], an intense effort was dedi- 
cated to investigate CoVs in bats, and many SARS-like 
viruses were detected and characterized in various species 
worldwide [4, 19]. Recent studies also provided evidence 
that bats may be the original ancestor hosts for SARS and 


identified in this study is shown in bold and assigned to the B-CoV 
genus. Bootstrap percentage values from 1000 resamplings are 
located at nodes of the tree 


MERS viruses, which recently caused pandemic and heavy 
consequences in humans [41-43]. Therefore, it is para- 
mount to investigate and characterize CoVs distribution 
and diversity in bat species, before their potential emer- 
gence in humans. 

Full-length genome analysis of the detected CoV should 
be attempted to thoroughly understand phylogenetic rela- 
tionships between SarBatCoV 1 and betacoronaviruses previ- 
ously described in European bats. Further studies including 
larger number of samples as well as continued surveillance 
on viruses circulating in Sardinian species are also needed 
to understand the role played by bats in the maintenance 
and transmission of beta CoVs. Monitoring and research 
on potential zoonoses are essential, and the prevalence 
and distribution of CoVs in Sardinian bats should be fur- 
ther investigated in urban and rural environment. However, 
without specific evidence, bats should not be considered as 
a major concern for public health, especially given the key 
ecosystem services they perform [44]. Bats play essential 
ecological roles in mosquito and pest control, as pollinators, 
and are keystone species in natural ecosystems. On top of 
this, bat species are often protected by national and interna- 
tional laws, as their conservation status may be vulnerable 
or threatened, and Sardinia was recently indicated as a major 
Mediterranean hotspot for bat biodiversity [26]. 
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g | KJ473819 BtMd-BetaCoV/HuB2013 


82 | 'KY770860 BatCoV isolate Jiyuan-84 Rhi_fer 


97 


DQ022305 Bat SARS coronavirus HKU3-1 


KT444582 Bat SARS-like CoV WIV16 
85 NC004718 SARS coronavirus 
97 'KC881005 Bat SARS-like CoV RSSHC014 


96 | 96 SarBatCoV1 


KC633199 Betacoronavirus BtCoV/Rhi_fer/It17/ITA/2009 
89 |- GU190215 Bat coronavirus R.blasii BM48-31/BGR/2008 
55 L GU190230 Bat coronavirus Rhi_ fer BM48-35/BGR/2008 
GQ404797 Bat coronavirus SLO1A0082/2008/SVN 


98 


KF500954 Betacoronavirus R.hipposideros/Italy/243585/2012 


Beta 


64 - KF500953 Betacoronavirus R.hipposideros/Italy/196814/201 1 
DQ249214 Bat coronavirus HKU4 isolate HKU4-1 
DQ249217 Bat coronavirus HKU5 isolate HKU5-1 
99 96 [- KF500941 Betacoronavirus N.noctula/Italy/206679-5/2010 
50 HQ184062 Betacoronavirus E.isa/M/Spain/2007 


75 


KC243390 Bat coronavirus BtCoV/8-724/Pip_pyg/ROU/2009 


JX869059 Human betacoronavirus 2c EMC/2012 2c 


70 


HQ184059 Betacoronavirus H.sav/J/Spain/2007 


de | KF500946 Betacoronavirus P.kuhlii/Italy/206645-53/2011 


100 


KF500950 Betacoronavirus P. kuhlii/Italy/330375-15/2012 


97 || KF500940 Betacoronavirus H.sav/Italy/206645-40/2011 
55 L KF500944 Betacoronavirus P.kuhlii/Italy/206645-29/201 1 
DQ249235 Bat coronavirus HKU2 isolate HKU2-1 
DQ249224 Bat coronavirus HKU6 isolate HKU6-1 


HQ184058 Alphacoronavirus P.kuh/Iprima/Spain/2007 
HQ184051 Alphacoronavirus N.las/C/Spain/2007 


Alpha 


60 KF500945 Alphacoronavirus P.kuhlii/Italy/206645-4 1/2010 
100 ' KF500945 Alphacoronavirus P. kuhlii/Italy/206645-41/2010 


0.2 


Fig. 2 Phylogenetic analysis based on 396 bp of the RdRp gene. Phy- 
logenetic tree of the partial RNA-dependent RNA polymerase gene 
(396 bp) of coronavirus strains found in bats and humans, showing 
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the assignment of the sequence detected in this study (SarBatCoV1) 
to clade 2b, within the B-CoV genus. Percentage values from 1000 
replicates of bootstrapping are located at nodes of the tree 
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